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Abstract 

We propose a new regularization method for constructing a shock wave type solution 
with nonsmooth front (interaction of shock waves) for quasilinear equations in the one- 
dimensional case. 



1 Introduction and Main Results. 

1. Consider the quasilinear first order equation 

L[u] =u t + (/(«))„ = 0, 



(1.1) 



where f(u) is a smooth function of at most polynomial growth, u — u(x, t), x £ M, t £ (0, T), with the initial 
value u(x,0) = u (x). The Hopf equation 



L H [u] = u t + {u 2 ) x = 0. 



(1.2) 



is a special case of this equation. 

It is known that equation ([0]) can have a discontinuous shock wave type solution even in the case of the 
smooth initial condition u°(x). Such piecewise smooth solutions are determined by integral identities, namely, 
the function u(x, t) is called a generalized solution of equation ( [L~l| ) in the domain Q C R if for all functions 
with compact supports <p(x, t) £ T>(SY) we have 



u(x,t)<p t (x,t) + f(u(x,t))<p x (x,t) 



dx dt = 0. 



(1.3) 



O. A. Oleinik |JJ proved the existence of a solution of the Cauchy problem for equation (LI) in the class of 
piecewise smooth functions provided that certain additional stability conditions hold. A theory of such problems 
for quasilinear equations and systems of equations is described in B - || 

However, the Cauchy problem can be posed for equation (LI) with the singular initial value u°(x) if its 
singularity is stronger than that of the Heaviside function 



1, t > o, 

0, £ < 0, 
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i.e. a jump. 

Thus, in Q - (9) the problem about the propagation of an infinitely narrow <5-soliton was considered. In 
the simplest case nonlinear waves of this type arise if we consider the weak asymptotics, as e — * +0, of the 
one-soliton solution 

wOM,e) = ych [—(x-vt)/e), 
of the Korteweg-de-Vries equation (KdV) 

LxdvVA =ut+ (u 2 ) x + £ 2 u xxx = 0. 

Up to Ox>'(£ 2 ), this asymptotics gives the infinitely narrow S-soliton 

u s (x, t) — AeS(x — vt), e — > +0, 

where S(x) is the Dirac 5-function, A = ^ / clT 2 (£) d£_ = 6^ and by O v >{e a ) we denote a distribution from 
D'(M.) such that for any test function ip(x) £ T> 

(Ov>(e a ),<p(x)) -0(e Q ), 

and 0(e a ) is understood in the ordinary sense. We stress that here and below we consider all distributions as 
distributions depending on the argument x. All other arguments are considered as parameters. 

Since u(x,t,e) — Ox>'(^) and s 2 u xxx = Ou'(e 3 ) as e — ► +0, the limit expression u e (x,t) was interpreted by 
V. P. Maslov and V. A. Tsupin fij, V. P. Maslov and G. A. Omelyanov Q as a generalized solution (asymptotic 
up to Oxj/(e 2 ) of the Hopf equation (1.2) which is the limit problem for the KdV equation. In 0, jij the 



corresponding generalized Hugoniot conditions similar to those at a shock wave front were obtained. Here the 
initial value u°(x) = AeS(x), e — > +0 is not a distribution but an asymptotic distribution (see M and below). In 



H asymptotic generalized solutions of the quasilinear equations (1.1) and systems of equations were considered 
in the form of infinitely narrow 5-solitons in the algebra of asymptotic distributions. 

In B the problem on propagation of infinitely narrow P-solitons was also considered. This is a new type of 



nonlinear waves which arise in solving the Cauchy problem for equation ( LI ) with the initial value in the form 
of the asymptotic distribution u°(x) — Uq(x) + g°(x)eP(x~ 1 ), e — > +0, where P(x~ 1 ) is the principal value of 
the function 

In |Io| - fl2|| problems on propagation and interaction of 5-w&ves for semilinear hyperbolic systems are 
studied, i.e. the Cauchy problem is solved for the initial value u (x) whose singularity is of the Dirac 5-function 
type. 

Generally speaking, the class of problems involving determination of singular solutions of quasilinear and 
semilinear equations can lead to the problem to define multiplication of distributions (generalized functions) 



and to construct associative algebras containing distributions. In particular, if we rewrite equation (1.1) as 
Ut + f'(u)u x — which is equivalent to the divergent form for the case of smooth solutions, then even for 
piecewise smooth solutions of this equation with jumps there arises a problem to define the product of the Dirac 
(5-function by the Heaviside function <5(^)_ff(^). 

These problems require a development of special analytic methods. There exist various approaches to their 
solutionpjllj, || -@. 

In |g7 and especially in || was developed a new analytical method, the weak asymptotic method, which 
enables investigation of the dynamics of propagation of various types of singularities of quasilinear differential 
equations and hyperbolic first order systems. The fundamental ideas of this method originate from the papers 
by Y. B. Livchak (T5J , Li Bang-He |l6), V. K. Ivanov JL7| (where the product of distributions was defined 
as the weak asymptotics of the product of approximations of distributions being multiplied with respect to 
the approximation parameter) and V. P. Maslov H (where the direct substitution of a singular ansatz into a 
quasilinear equation was used). 

In the present paper the weak asymptotics method is applied to the investigation of the dynamics of both 
propagation and interaction of initial discontinuities, i.e. shock waves. Further we describe the essence of the 



weak asymptotics me thod and give the definitions 1.2, 1.3 of a generalized asymptotic solution and a generalized 
solution to equ ation ([L~l|). 



In Theorem |l . l| and Corollary L3, Corollary L4, using the weak asymptotics method, we wri te out a system 



of equations to determine the dynamics of propagation of one shock wave for equation (1.1) and the Hopf 



equation (b2), i.e. in the class of piecewise smooth functions we solve the Cauchy problem with the following 
initial value: 

U {X '-\ <(x), x > xo, (L4) 
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where xq is the initial position of the shock wave front, Uq(x), e (x) > are smooth functions. Analogous 
results for strictly hyperbolic first order systems were derived in |J. The proof of Theorem 1.1 is given in 
Section ||. 



In Theorem 1.3 we write a system of equations to determine the dynamics of propagation and interaction 



of two initial discontinuities (shock waves with constant amplitudes) for the Hopf equation (1.2) is considered 



since it is a typical model of a nonlinear hyperbolic equation. In Theorem 1.5 the results obtained are extended 



to the case of equation (1.1). Namely, we solve, in the sense of Definition 1.3, the Cauchy problem for equations 



(1.2) and (1.1), in the class of piecewise constant functions, with the following initial condition: 



M 



■ ei 
u 



■ e 2 , 
- e 2 , 
"o, 



x ^ x^ , 

x\ < X < X 

X > Xn 



2' 



(1.5) 



where u$, e±, e 2 are constants. The solution of the Cauchy problem is given by the formula which can be used 
for all t > 0. 

To prove Theorems 1.3, 1.5 we construct asymptotic solutions of equations (1.2) and ( |1 . 1[) in the sense of 
Definition 1.2. These asy mptotic solutions are given by Theorems 1.2, 1.4 proved in Sections |4j, ^. 
In Theorems 1.2 , 1.4 we describe the dynamics of the shock wave merging process. 

In what follows we assume that f"(u) > and the solutions u(x, t) are piecewise smooth functions satisfying 
the Oleinik stability conditions at every point (x, t) of a discontinuity line 



u(x- 0,t) > u(x + 0,t). 



(1.6) 



The o ne-p hase and two-phase generalized solutions constructed in Theorems 1.1 , 1.3 



L5| in the sense of 

Definition L3 are generalized solutions in the sense of the standard definition by the integral identity (|l.4| ). 
Thus, provided that the stability conditions are fulfilled, the Oleinik uniqueness theorem holds for them M. 

Note that our methods are also applicable for solving a similar problem in the case of systems of quasilinear 
first order equations. 

2. The weak asymptotics method. When solving the problem on propagation and interaction of singular- 
ities of quasilinear equations one has to extend the Schwartz distribution space. Therefore, in || we constructed 
an associative algebra of asymptotic distributions with the identity and free of zero divisors generated by the 
linear span of (one-dimensional) associated homogeneous distributions. The elements f£(x) from the algebra of 
asymptotic distributions are defined as weak asymptotic expansions of elements f*(x,e) from the linear span 
of approximations of one-dimensional associated homogeneous distributions as the approximation parameter e 
tends to zero. Each element of this algebra has a unique representation in the form of the asymptotics (in the 
weak sense) whose coefficients are distributions. 

The product of associated homogeneous distributions is defined as a weak asymptotic expansion of the 
product of approximations of multiplied distributions as e — > +0. This product is an element of the associative 
algebra of asymptotic distributions. 

In particular, infinitely narrow 5 and P-solitons are elements of the algebra of asymptotic distributions. 

Now turn to the description of our technique omitting the algebraic aspects which are given in detail in 0j 
and 111. 



To study the interaction of (two) singularities of equation (IT), we solve the initial value problem with 
distribution initial data 



i*°(x) = u° (x) + J2e Q J {x)s 1 (x - Xl ), 



where Uq(x), e|? (x) are smooth functions, Sj(£) are distributions (generalized functions) or asymptotic distri- 
butions (see ||, jl!|) and Xj are constants. 

The solution of this initial value problem is found in the form of the singular ansatz 



u*(x, t) — u (x, t) + 



j(x,t)Sj(x - <j)j(t)) 



(1.7) 



where uo(x,t), ej(x,t), (f)j(t) are functions to be found and in the general case (f>j(t), ej(x,t) and Sj(£) may 
depend on the small parameter e. The singular ansatz u*(x,t) belongs to the associative and commutative 
differential algebra of asymptotic distributions, However, actually, we do not use this fact in this paper directly, 
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since to define the associative and commutative product of distributions, we need to calculate all terms of the 
asymptotic expansion (i.e., to calculate the weak asymptotics of the product of approximations up to Op'(e°°)). 
But in view of the above remarks and Definition 1.3 below, to construct a generalized asymptotic solution, it is 
sufficient to calculate the weak asymptotics of the above-mentioned product of approximations up to Ox>/(e°°), 
N = 1 or N = 2 (see below). 

To study the interaction of two nonlinear waves we assume: 1) Sj(^) = j = 1,2 in the case of the 

5-waves; 2) Sj(£) = H(£), j = 1,2 in the case of the shock waves; 3) Sj(£) = £<$(£), j — 1,2 and Sj(£) = eif (£), 
j = 3, 4 in the case of infinitely narrow 5-solitons; 4) Sj(£) — eP(^ _1 ), j = 1, 2 in the case of infinitely narrow 
P-solitons. 

To find the solution of the form u s (x, t) we construct the smooth ansatz 



*(x, t, s) — Uq(x, t) + ej(x, t)sj (x — (f>j(t, s),s), 



(1.8) 



which is a smooth approximation of the singular ansatz ( |l.7| ). Here Sj(£, e) is a smooth approximation of the 
distribution or asymptotic distribution Sj(£) and e is the approximation parameter. 

To construct the approximations Sj(£,e) we use the fact that a correspondence can be set up between each 
distribution (generalized function) fix) € V and its approximation: 



f(x, e) = f{x) * K{x, e) = (fit), K(x -t,e)), e> 0, 



(1.9) 



where * is a convolution, the kernel K{x,e) — -w(-) is a <5-type function such that < ui(z) £ C°°(M.), 
J u>(z) dz = 1, u>(z) has a compact support or decreases sufficiently rapidly as \z\ — > oo, for example, |w(z)| < 
(7(1 + \z\)~ N for all positive integers N @ ch.I, §4.6.]. 
For all test functions <p(x) € V we have: 

Urn {f(x,e),tp(x)) = (f(x),ip(x)). 

s— >+0 



For example, for an approximation of (^-function we have from (1.9): 

6{x,e) = -uj(-). 
e e 



(1.10) 



For an approximation of the Heaviside function H(x) we have from (1.9): 

H(x, e) = H(x) * -«(-) = / w(- - t) dt. 

£ £ Jo e 



Hence we find 



H(x,e) = ojo(-) 



u(ri) di], 



(1.11) 



where wq(z) € C°°(R), lim z ^ +00 cJo(^) = 1, hm 2 



,u (z) = 0. 



The smooth ansatz u(x,t, e) is substituted into the equation (1.1) and the weak asymptotics (in terms of 
the space of the distributions V) up to 0-p>ie N ), where N = 1 for delta and shock waves, N — 2 for infinitely 
narrow solitons, is found of the left- han d side of this equation as e — > +0. We define the generalized solution 
(asymptotics solution) to equation (IT) as the weak asymptotics u B (x,t) of the smooth ansatz u(x,t, e), as 
e-> +0. 

The key role in solving the problem of the interaction of nonlinear waves of the equa tion (1.1) plays the 
construction of the weak asymptotics of the smooth function from the smooth ansatz (1.8), which can be 
represented in the form 



f(u (x,t) +y]ej(x,t)sj(x - <f>j(t, e),e) 



3 = 1 
2 



= f(uo(x,t)) + ^ Bj (up (x,t), e x (x,t), e 2 (x,t), p) Sj (x - <j>j(t)) + Ov(e N ), 



(1.12) 
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where the estimate Oj)/(e ) is uniform with respect to tl>(t, e), where ip(t, e) = fait, s) — </>i(i, e), p 



V-(M) 



, and 



the interaction switches Bj(uo(x,t),ei(x,t),e2(x,t),p), j = 1,2 which are smooth functions, can be computed 
explicitly. 

Expansion (1.12) allows one to separate the singularities of the ansatz. Therefore, substituting expansion 



(1.12) and the singular ansatz (1.7) into a quasilinear equation and setting equal to zero the coefficients of the 
different powers of the small parameter s and of the linear independent distributions we obtain a system of 
equations (in particular, the Rankine-Hugoniot type condition) which describes the dynamics of singularities 
and defines the smooth functions uq (x, t), ej(x,t), 4>j{t), j — 1,2. 

If the prop agat io n of one singularity is studied we set e^x) = in the initial value u*°(x), and we set 
e 2 (x,t) = 0in (p|), (pj|). 

3. The weak asymptotics method and shock wave s. It fol lows from (1/7), (|l.8| ) that to study 
propagation of a shock wave, i.e. to solve the Cauchy problem (1.1), (1-4), we must seek a solution in the form 
of the singular ansatz 

u* e {x,t) =u (x,t) + e(x,t)H(x -</>(*)) (1-13) 



and substitute into equation (1.1) the smooth ansatz 

u*(x, t, e) = uq(x, t) + e{x, t)H{—x + </>(i), e), 



(1.14) 



where Uo(x,t), e(x,t), <fi(t) are functions to be found, (f)(0) = Xq and is the approximation (1.11) of 

the Heaviside function H(^). 

It should be noted that an the algebraic approach to the construction of solutions of the shock wave type, 
that is, the method of direct substitution the singular ansatz (1.13) into a quasilinear equation in divergent 
form was first used by V. P. Maslov ||. His approach was based on the fact that a smooth function / of the 
sum uq(x) + e(x)H(x), where uq(x), e{x) are smooth functions, can be represented in the form as its argument: 

f(u (x) + e{x)H{x)) = f(u (x)) + [f(u (x) + e(x)) - f(u Q (x))]ll(x). 

After substituting this expression into a quasilinear first order equation written in the divergent form, and 
differentiating, one can obtain a linear combination of the unit, the Dirac (5-function and the Heaviside function 
with smooth coefficients. Setting this expression to zero and separating the singularities, we obtain the relations 
which determine the solution. 

In order to investigate the interaction of two shock waves, that is, to solve the Cauchy problem (1.1), (|l.5|), 
we should seek a solution in the form of the singular ansatz 



u* e (x,t) = u (x,t) + ^ \ek{x,t)H( - x + (p k (t,e)) 

k=l 

and substitute the smooth ansatz 

2 

u*(x,t,e) = u (x,t) + \e k (x,t)H k (- x + cj) k (t,e),e) 

k=l 



(1.15) 



(1.16) 



into equation (1.1). Here u (x,t), ek(x,t) are smooth functions to be found, lim e ^ + o 0fe(O, e) = x { ^, iffc(£,£j 
are approximations ( 1.11 ) of the Heaviside function H(£), k — 1,2. 

Definition 1.1 Let f(u) be a smooth function, and let u*(x,t) is an asymptotic distribution (singular ansatz) 
(1.13) or (1.15), and it*(x, i, e) is a smooth ansatz (1.14) or ( 1.1€| ). 

By the Ox>< (e)-substitution we call the weak asymptotics of the expression / (u* (x,t,e)) as e — ► +0 computed 
up to Ox>i{e). We denote this substitution f(u*(x,t)). 

Now we present the main lemma which gives an asymptotic expansion of the type ( 1.12| ) for the case of 
shock waves. The proof of Lemma 1.1 is given in Section 0. 



Lemma 1.1 Let 



H(-x,s) 



wi(r))dr], 



H(—x + a, e) 



w 2 (f?) dr\ 
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be approximations of the Heaviside functions H(—x), H(—x + a), respectively, the approximating functions 
uJk(z) £ C°°(]R) are nonnegative and either have compact supports or decrease sufficiently rapidly, as \z\ — ► oo 
and J tOk{z) dz — 1, k = 1, 2. 

Let f(u) be a smooth function of at most polynomial growth and let uo(x,t), e±(x,t), e 2 {x,t) be bounded 
functions. Then we have the asymptotics: 

f(u (x,t) + e 1 (x,t)H 1 (-x,e) + e 2 (x, t)H 2 (-x + a, e)) = f(u (x,t)) 
f(u Q (x,t)+e 1 (x,t))~f(u (x,t)) H(-x)+ f(u (x,t)+ex(x,t))- f(u (x,t)) H(~x + a) 

-0, 



+ Bi(x, t, -)H{-x) + B 2 (x, t,--)H(-x + a)+ O v ,(e), e - 

where the estimate Ot>>(s) is uniform with respect to a. 

The functions Bk(x,t, p), k = 1,2 called interaction switch functions have the following form 

Bi(x,t,p) = / f'(uo{x,t) + ex (x, t)uoi (—»?) + e 2 (x,t)uo2(-r) + p)) 

-f'(u (x,t) + ei(x,t)u i(-r))) e 1 (x,t)u 1 (-i]) dn, 
B 2 (x,t,-p) = J f'(u (x,t) + ei(x,t)cJoi(-T) - p) + £2 (x, t)uj 02 (-??)) 

-f'(uo(x, t) + e 2 {x, t)u3{H2{—rj)) e 2 (x, t)uj 2 (-Tj) dr\. 

In addition, the interaction switch functions satisfy the following relations: for each p £ M. 

B 1 (x,t,p) + B 2 (x,t, -p) = f(u (x,t) + ex(x,t) + e 2 (x,t)) - f(u (x,t) + e 1 (x,t)) 

- f(u (x, t) + e 2 (x, t)) + f(u (x, t)), 
lirn /9 _ >+0o Bk(x,t,p) = f(u (x,t) + ei(x,t) + e 2 (x,t)) - f(u (x,t) + ei(x,t)) 

- f(u (x, t) + e 2 (x, t)) + f(u (x, t)), 

limp^-oo Bk(x,t,p) = 0, fc = l,2. 



(1.17) 



(1.18) 



(1.19) 



Remark 1.1 For simplicity, we have selected the approximating functions u>k{z) such that they either have 



compact supports or decrease sufficiently rapidly as \z\ — > oo, for example, |wfc(z)| < CV-(1 + \z\) 
positive integers N, k = 1, 2. Therefore, we have 

u<3k(z) = f* o0 UJk(v)dri = l + 0(z- N ), z^+oo, 
u> 0k {z) = 0(\z\- N ), z->-oc. 

Consequently, applying the Lagrange theorem we obtain the estimate 

[f'(u + eiu>oi(-r)) + e 2 u a2 (-r/ + p)) - f'(u + ei^oi (-»?)) e.\Ux{-rj) 
= /"(K + eiW i(-?7)] + ®e 2 uj 02 (-r) + p))eie 2 ui(-r))uj 02 (-r/ + p), < 6 < 1. 



-N 



for all 



It follows from this estimate and (1.18) that 

Bx(x,t,p) = f(u Q + e x + e 2 ) - f(u + ei) - f(u Q + e 2 ) + f(u ) + 0{p~ N ), p -> +oo, 

B x (x,t,p)=0(\p\- N ), p^-oo, 
where N = 1, 2, . . .; if po is a constant then 

Bx(x,t,p) = B 1 {p ) + O{p- p ), p^po- 



Corollary 1.1 Let f{u) be a smooth function of at most polynomial growth and let uq{x, t), e(x, t) be bounded 
functions and H(~x, e) be approximations of the Heaviside functions H(—x). Then 



f(u (x, t) + e{x, t)H(-x + 4>{t),e)) = f(u (x, t)) 

+ \f(u (x,t)+e(x,t))-f(u (x,t))]H(-x + <t>(t)) + O v ,(£), e 



(1.20) 



G 



If we set e\{x,t) = e(x,t), e 2 {x,t) = in (1.18) then B\{x,t,p) — B 2 {x, t, — p) = and expansion (1.2C) 
immediately follows from (1.17). 



Remark 1.2 In view of the obvious relation 

[H(x,s)] r = H(x)+O v >(e), e -> +0, r = 2,3,. 
for the case when f(u) is a polynomial, formula (1.2C) was proved in j^, 8.1.]. 



(1.21) 



Corollary 1.2 Let 



H(-x,s) 



w\{t})dr], H(—x + a, e) 



be the approximations of the Heaviside functions H(—x), H(—x + a), respectively, where the functions u> k {z) £ 
C°°(]R) either have compact supports or decrease sufficiently rapidly for \z\ — > oo, and J Wk(z) dz = 1, k = 1, 2. 
Then 



H(-x, e)H{-x + o,e) = B x {~)H(-x) + B 2 (- -)H{-x + a) + Op/(e), e -» +0, 

where the estimate Op/(e) is uniform with respect to a. 

The functions B k (p) called interaction switch functions, have the form 

B i{p) = f-oc (&i * w 2 )(r/) efy, 
#2 (~p) = JT£ (wi * w 2 )(r/) efy, 

where &{rj) — w{—rj), * is ifte operation of convolution. 
In this case, 

B 1 {p) + B 2 {- P ) = l, peM, 
B k (oo) = lira B k {p) = 1, 

B fe (-oo)= lim B fc (p) = 0. 

p — ► — OO 



(1.22) 



(1.23) 



(1.24) 



To prove Corollary L2| we find the asymptotics of the expression 
H(—x, e)H{—x + a, e) 



qquad 



= ^(f(H 1 (-x,£)+H 2 (-x + a,e))-f(H 1 (-x,e))-f(H 2 {-x + a,e)) 



where f(u) = u 2 



Setting in (1.17) Uo(x,t) = 0, e±(x,t) — e 2 (x,t) — 1 and taking into account relation (1.21), we obtain 



( |1.22p . In this case it follows from ( |1.18| ) and ( |l.ll| ): 

B k ({-lf- l p) = I c 3- fc (-r/+(-l) fc -V))^(-^^ 

LOz-k{v) dv ujk(-ri) dr), k = 1,2. 

After the substitution v — z — £ in the second integral, we reduce B k (fl) to the form 

,(-i) fc -V 



■ fc-i 



Wfc(-0 w 3-fe(2 - d£ 



dz 



and thus obtain (1.23) 



In this case (1.19) implies (1-24). 

4. Generalized solu ti ons. Let us define a ge neralized discon tinu ous solution and rules for substituting 
ansatzs of the type ( 1.13 ), ( 1.15| ) into equation (1.1) (see Definition 1.1). 
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Definition 1.2 Let u*(x,t, e) be the smooth ansatz (1.14), ( 1.16| ). We call the corresponding asymptotic 
distribution (singular ansatz) (1.13), (1.15), u*(x,t), by a generalized asymptotic shock wave type solution of 



the equation L[u] = 0, for t £ [0, T], with the initial condition u®*(x), if the following relations hold 



L[u*(x, t, e)] 
u* E (x,0) 



Ov>(e), 
u°*(x) 



+ O v ,{e), 



where the first estimate is uniform with respect to t 6 [0, T\. 



Definition 1.3 Let u*(x,t) be an asymptotic of the order of Ou'(e) solution of the equation L[u] — with the 
initial condition vP e *{x). By a generalized solution of this equation with the initial condition vP e *(x) we call the 
weak limit 

u*(x,t) = lim u*(x,t). 

£—>+0 

After substitution of expansion (1.20) and the singular ansatz (1.13) (or expansion ( 1.17| ) and the singular 
ansatz (1.15)) in equation ( [l.l| ), to separate the singularities of the left hand side of equation ( [l.l| ) and to 
obtain a system of equations which describes the dynamics of singularities and determines the smooth functions 
uo(x,t), efc(cc,i), 4>k(t), k — 1, 2 we use the following Lemma. 

Lemma 1.2 If A(x), Bk{x), Ck{x) are smooth functions k — 1,2 and a > 0. Then 

A{x) + Bx{x)0{x) + C x {x)5{x) + B 2 {x)9{x - a) + C 2 {x)5{x - a) = 

if and only if 



A(x) = 0, 

A(x)+B 1 (x) = 0, 

A{x) + B x {x) + B 2 {x) = 0, 

Ci(0) = o, 

C 2 (a) = 0. 



for x < 0, 

/or0 < x < a, 
for x > a, 



5. Propagation of shock wave of equation ( |l . l| ) . 

Let us consider the propagation of the single shock wave of the equation L[u] = u t + (f(u)) x = with the 
initial value ( |1.4| ) . In the framework of our approach we shall seek the solution in the form of the singular ansatz 
( 



1.13) which coincides with the initial value u*°(x) for t — 0. This means that one substitute the smooth ansatz 

1.14) approximating ( 1.13| ) into the equation L[u] — and find the weak asymptotics, up to OV(e), of the left 
hand side of the equation, as e — > +0. 



Theorem 1.1 Consider the Cauchy problem (1.1), (1.4), where f(u) is a smooth function, f"(u) > and 
Uq(x), e°(x) > are smooth functions. Suppose there exists T > 0, such that 



inf 

€<0(o) 



- inf >°WV(6)] < T-\ 
A0+eV)l7>°(0+e°(£))] < T-\ 



(1.25) 



Then in the sense of Definitions 1.2, 1.3 equation (1.1) for t G [0, T] has a discontinuous solution of the 
form ( 1.13 ) 

u*(x, t) = u (x, t) + e(x, t)H(-x + 4>(t)), 
if and only if the unknown functions Uo(x, t), e{x, t), <p(t) G C°° satisfy the following system of equations 

(1.26) 



L[u (x,t)] = 0, 

L[uq(x, t) + e(x, t)} = 0, 

d(j)(t) _ f (up (z,t)+e(x,t)) -/ («o(g,t)) 

dt ~ ^) 



x > <t>(t), 

X < (f)(t), 



where u*(x,0) = [uo(x.t) + e(x,t)] 
x > 0(0), (f>(0) =x . 



u^(x) + e°(x) for x < 0(0) and u*(x,0) = uo(x,t) 



Uq(x) for 



8 



Corollary 1.3 In the case of the Hopf equation Lh[u] = u t + (u 2 ) x — system ( 1.26| ) has the fo 



L H [u (x,t)} = 0, 
L H [u (x,t) +e(x,t)] = 0, 

^ = 2u a (cf>(t) 1 t) + e^(t),t). 



x > 0(t), 

x < <j>(t), 



(1.27) 



The last equations in systems (1.26) and (1.27) represent the Rankine-Hugoniot conditions along the dis- 
continuity line x — ip(t) for shock waves || Ch. 4, §1.2.]: 



dx 
~db 



[/(«)] 



(1.28) 



where 



H = u(tp{t)+0,t)-u(<p{t)-0,t), 
[/(«)] = /(t*(¥»(t)+0,t))-/(u(^(t)-0,t)) 



are jumps of functions m(x, i) and /(u), respectively, along the discontinuity line x = ip(t). In our case, u = uq + c 
to the left of the discontinuity, and u — uq to the right of the discontinuity. 



Corollary 1.4 Solution (1.13) constructed in the sense of Definition 1.2 to the Cauchy problem for equation 
(1.1) with the initial value (1.4), that is, a function of the form u*(x, t) = uo(x, t) + e(x, t)H(~x + 4>(t)), whose 
smooth components satisfy system of equations (1.26) from Theorem 1.1, is the unique generalized solution of 



the same Cauchy problem in the sense of the standard definition by integral identity fll.q ). 

Proof. To prove that the solution of the Cauchy problem in the sense of Definition |l.3| is the solution in 
the sense of the standard defin ition by the int egral identity (1.3) we remind that in view of Theorem 1.1, the 
approximation u*(x, t, e) ( 1.14 ) of the solution ( 1.13 ) satisfies the equation L\u*(x, t, e)] =0-D'(e). Let us apply 
the left-hand and right-hand sides of this relation to an arbitrary test function (p(x,t) £ £1 C R 2 . 

Since for e > the functions u*(x, t, e) is smooth, integrating by parts, we obtain 

u*(x,t,s)cp t {x,t) + f(u*(x,t,s))ip x (x,t) dxdt = (Ov(e),ip(x,t)) = 0(e). 

Since the limit u*(x, t) of the family u*(x, t, e) as e — > +0 in the V sense is a bounded locally integrable function, 
we have J n 
identity ( |l~3 



-0, which coincides with the usual integral 



u*(x, t)<pt{x, t) + f(u*(x, t))ip x (x, t) dxdt = 0, as e 
) (see fj- |). 

Due to our assumptions, the stability condition ( |l.6| ) is satisfied and by Theorem of O. A. Oleinik this 
solution is unique 0. 

6. Interaction of the shock waves of the Hopf equation (1.2). 

To study the interaction of two shock wave s of th e Hopf equation (h2) Lh [u] = u t + (u 2 ) x = we shall seek 
its solution in the form of the singular ansatz ( 1.15| ) which, for t = c oincides with distribution ( |l.5[) . 

First we construct an asymptotic solution of problem (1.2), (1.5) in the sense of Definition |l.2| , which is 
given by Theorem 1.2 proved in Sections H. 



Theorem 1.2 For all t S [0, +oo) in the sense of Definition 1.2 there exists the asymptotic solution 

2 



u*(x, t, e) =u +} j e k H k ( - x+ 4> k {t, e), e) , 



fc=i 



for the Hopf equation (|l.2| ) such that: 

1) for allt > 0: L H [u*(x,t,e)} = Ov{e), 

2) u*°(x) - u*(x,0,s) = O-D'ie), where 



u*°(x) = 



- ei + e 2 , 
uq + e 2 , 
""0, 



x < ar", 



x > x 



'■2- 



x\ < x® and uq, e\ > 0, e 2 > are constants; 
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3) for t e (0, t*) 



lim u*(x,t,s) = uq + e k H( -x + 4>ko{t)), 

e— >+0 



k=l 



where 



fco (O) = x° k , k = l,2, t* 
4) for t G (**, oo) 



^io(i) = lim 0!(t, e) = 1O (O) + (2u + e x + 2e 2 )t, 
<ho{t)= Urn 6 2 (t,e) = (f> 20 (0) + {2u + e 2 )t, 

no 

1 is the shock wave merging time; 



ei+e 2 



lim u*(x, t, e) = u + (e\ + e 2 )H( — x + </>_(£)) , 

e— >+0 x ' 



wht 



<j>-{t)= lim <t> k (t,e) = x* + {2u (s + e 1 + e 2 ){t-t*), fc=l,2, 



2;* = 0fco(**) 



(2u + ei + 2c 2 )x\ - (2u + e 2 )x\ 



Here 



H k (-x,e) 



ei + e 2 



u>k(v)dr), 



are approximations of the Heaviside function H(—x), where tOk(z) are smooth functions which either have 
compact supports or decrease sufficiently rapidly as \z\ — > oo, J u>k(z) dz — 1, k = 1, 2. 
The phases 4>k{t,e) have the form 



r _Mj ' 



0jfe(f,e) = <feo(*) +VoW^*i(t) 
where ipo{t) — (f> 2 o(t) — 4>w{t) and the perturbations of phase functions are given by formulae ( 4.61 ): 

0*i(r) - (-l)^ 1 , 2 ;"\ f (l - Bi(p(t0)) 



According to (1.23), 



Bi{p) 



(fii *u 2 )(r])dr), 



and p = p(t) = [<f> 2 (t, e) — 4>i(t, e)]/e is a solution of the differential equation with the boundary condition ( 4. 52] ), 

= 2 Bl (p)-l, 



Note that, in fact, Theorem 1.2 is based on the remark that the variable t is considered as a parameter. If 
we considered the variable t as a variable having the same rights as the variable x, we also had to calculate a 
weak asymptotics with respect to the variable t. In this case we would obtain Heaviside functions instead of 
interaction switch functions and could not obtain a formula that were uniform in t for the solution u*(x, t, e). 

Theorem 1.2 implies the following result. 



Theorem 1.3 The Cauchy problem for the Hopf equation ( |l.2| ), (1.5) in the sense of Definition 1.3 has a 
discontinuous solution of the form ( 1.15| ) 



u*(x,t) = u a + ^ [efc-H"( - x + 4>k{t)) , 
fc=i 

where 

<M*) = 0io(O) + (2«o + ei + 2e a )t-e 2 (t-t*)fl"(t-f), 
Mt) = (j) 2Q {0) + {2u Q + e 2 )t + e 1 (t-t*)H(t-t*) 

are phases. 

This solution of the initial value problem ( |l.2[ ), (|l.5| ) is the unique generalized solution of this Cauchy problem 
in the sense of the integral identity (1.3). 



(1.29) 
(1.30) 
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Proof. It follows from Theorem |1.2| that distribution (1.2S) is a generalized solution of the initial value 
problem ( |l.2[ ), (1.5) in the sense of Definition [O . 

In view of Theorem 1.2, the approximation of the solution (1.29) uniformly in t satisfies the equation 

L[u*(x, t, e)] — Ox)'(e)- By applying the left- and right-hand sides of the last relation to an arbitrary function 
<p(x,t) £ we obtain 



u*(x, i, e)(ft(x, t) + f(u*{x, t, ey\ip x {x, t) dxdt = O(e). 

By integrating by parts for e > and then passing t o the limit as e — > +0, just as in the proof of Theorem |Q] , 
we prove that the functions u*(x,t) determined by ( 1.29| ) satisfies the integral identity 



u*(x,t)tp t (x,t) + f(u*(x,t))ip x (x,t) 



dxdt = 0. 



Since for the piecewise constant solution u*(x,t) the stability condition > 0, k = 1,2, holds, this solution is 
unique by the Theorem of O. A. OleinikJpJ . Theorem 1.3 is proved. 

The solution described in Theorem 1.3 is of rather simple natural form. In fact, it is sewed together of the 
following parts: two shock waves till the time instant t* and one shock wave for t > t* . We cannot say anything 
about the equation at the point x = x* , t = t* , but we see that the trajectories of discontinuities at this point 
are continuous. Moreover, one can easily see that the function ( 1.30| ) itself is continuous w.r.t. to t in the weak 
sense. Thus we have the following "naive" situation: everywhere outside the curves x = <f>k(t), k = 1,2, the 
solution (1.30) ca n be f ound as a function satisfying the equation L[u] — 0, while the curves x = 4>k(t), k = 1, 2, 
and the solution (|l.30|) are determined by the Rankine-Hugoniot conditions for t ^ t* and the condition that 
they are continuous at the point t = t* . 

However, as we shall see now, one can construct a great deal of such points. As an example, we consider the 
following Cauchy problem for the Hopf equation: 

u t + {u 2 ) x = 0, u(x,t)\ t=Q = H(-x - 1) - H(x - 1). 

One can easily verify that for t < t* = 1 the solution u(x, t) consists of two shock wave approaching each other: 

u(x, t) = H(-x + 1 - 1) - H(x + 1 - 1). 



By using the techniques used in Theorem |1.2| for calculating the solution u*(x,t), we see that the solution 
satisfying the integral identity for t > t* = 1 (this solution is one of "naive" solutions) is the stationary solution 
of the Hopf equation 

u(x,t) = H(-x) - H(x) 

consisting of two shock waves sticking together. 

Another different "naive" solution can be constructed as follows. Let us consider the solution of the Hopf 
equation with the initial condition 



W(x,t) 



t=o 




The solution W(x,t) of this problem has the form 

W(x, t) = H(-x - <j){t)) - H(x - <p{t)) + 



1 + 2at 



x < -1, 
-1 <x < I, 
x > 1. 



[I -H{-x-<f>{t))-H{x -<!>{€))]■ 



Its plot consists of two horizontal lines (y = 1 to the left of the point x_ = — <f> and y = 1 to the right of the 
point x+ — 4>) and the inclined line y = ax/(l + 2at) between the points a;_ and x + . 

At the points x± the solution has a jump whose coordinates x± — ±4>(t) can be calculated from the 
Rankine-Hugoniot conditions 



(-l) 2 -(aa;/(l + 2 a t)) 2 



-1 - ax/(l + 2at) 



x—(p 



l + 2at 



1, 



0(0) = 1. 



Assume that a > 1. Then, obviously, 4>t > on some interval [0, to]; to > 0. 
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We consider the limit of the function W(x,t) in T>'(M. X ) as t 
have 

lim (W(x,t),(p) — lim 

t->-i/(2o)+o t-»-i/(2o)+0 Ll + 2at 

By using the L'Hospital rule, we obtain 



-l/(2o) + 0. For any test function ^?(x) we 



z^(x) dx + 



dec 



<£>(x) dx 



lim (W, <p(x)) 

t^-l/(2a)+0 



lim 

f->-l/(2a)+0 



2 a 



dec 



y(cc) dec 



Consider the function 



U(x,t) 



(H(-x) - H(x))cp(x) dx. 



u(x,t) 0<t<t* 
W{x,t-l/{2a)-t*) t>t*. 



x± 



This function is continuous in T>' and satisfies the classical Hopf equation everywhere except for the discontinuity 
trajectories 

±(t- 1) t<t* = 1, 

±0(x,<-l/(2a)-t*) i>l, 

which arc continuous curves. Thus this is the second "naive" solution of the problem of discontinuity merging. 
In this solution, the discontinuities seem to interact elastically and then to run away after the interaction. 

Thus Theorem [Tl^ is not useless. This theorem allows one to choose a "true naive" solution among all 
possible solutions. 

The example of solution W(x, t) is of independent interest, since this is one of few solutions with nonmonotone 
motion of its singularity []. Namely, if a > 1, then, as was already pointed out, the discontinuities run away from 
each other but, starting from the time instant at which the inclined line y — ax /(I + 2at) enters the interior of 
the strip y = ±1, the directions of their motion are changed to the opposite, and the discontinuities merge into 
the stationary solution 

u(x,t) = H(-x) - H(x). 

Moreover, at the time instant of merging, we again can "paste in" a W(x,t) type solution at the discontinuity 
point, repeat this process, and obtain the picture of i-periodic elastic interaction of shock waves. 
7. Interaction of shock waves of equation ( |l . lj ) . 

Study the process of interaction of two shock waves of the equation u t + {f(u)) x — 0, where f(u) is a convex 
(downwards, i.e., f"(u) > on the range of the solution u) smooth function. To this end, consider the Cauchy 
problem with the initial value (1.5). As in the case of the Hopf equation, substitute the smooth ansatz (1.16), 
where </>i(0) = x®, ^2(0) = x^ and e\, > 0, approximating the singular ansatz ( 1.15| ), into equation (1.1) 
and find the weak asymptotics of its left hand side. 

Eventually, we construct the asymptotic solution of problem (p..lj), (1.5) in the sense of Definition |1 .2] , which 
is given by Theorem 1.4 proved in Sections ||: 



Theorem 1.4 For all t G [0, +00) in the sense of the Definition |1.2j there exists the asymptotic solution 



u*(x, t, e) = u + 2J e k H k ( - x + (f> k (t, e), e) , 



k=l 



of the equation (1.1) such that: 

1) for all t > 0: L[u*(x,t,e)] = Ov(e), 

2) u*°(x) - u*(x,0,e) = Ot>'(e), where 



u*°(x) 



x i < x 2 an d uq, ei > 0, e2 > are constants; 




x < X®, 

X-y X ^ 1 

x > cc!], 



1 G. A. Omel'yanov helped the authors to clear up this situation. 
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3) for t e (0, t*) 



lim u*(x,t,s) = uq + e k H( -x + 0fcoft)), 

e— >+0 



k=l 



wht 



Mt)= lim Mt,e) = ^ 10 (0)+ /(no+ei+e : i ) " /(Mo+e2) ^ 



(0)=^„ fc=l,2, 



t = eie2 



'20 



(0) 



/(M +e 2 )-/(Mo)^ 



r _ ~0 
L 2 x l 



e2/("o + ei + e 2 ) - (ei + e 2 )/(u + e 2 ) + eif(u ) 



is the shock wave merging time; 
4) for £ G ft*, oo) 



where 



lim u*(x, £, e) = u + (ei + e 2 )H( — x + (j)-(t)) , 

e— >+0 



Lft) = lim ^(t,e) = + /( "° + 61 + 62) " /(U0) ft ~ f), * = 1,2, 
e^+o ei + e 2 



/(u + ei + e 2 ) - /(ito + e 2 ) 



e 2 a;2 



/(u + e 2 ) - f(u ) 



Here 



e 2 /(Mo + ei + e 2 ) - (ei + e 2 )/(w + e 2 ) + eif(u ) 
H k (-x,e) = / d?7, 



are approximations of the Heaviside function H{—x), where u>k(z) are smooth functions which either have 
compact supports or decrease sufficiently rapidly as \z\ — > oo, J tOk(z) dz — 1, k = 1, 2. 
The phases (f>k(t, e) are set in the form 



(j>k(t,e) = <f>ko(t) + -0oft)<?W(i~) 



*n(t) ' 



where ipo ft) — </>20 ft) — </>io ft) an d the perturbations of phase functions are given by the formulas 
Mr) 

T 

_ , ^k-i e 3-fe / /("o + ei + e 2) - /(mq + ei) - f(u + e 2 ) + f(u ) - £i(p(t')) ^ 
T /("o + ei + e 2 ) - /(w + e 2 ) e 2 - /(w + e 2 ) - /(u ) ei 

Bi(p) = / /'(u + eiw i (-??) + e 2 cj 2(-^ + p)) - /'(u + ei^oi (-77)) eiwi(-r]) drj, 



where p = p{r) has the same meaning as in Theorem 1.2 and satisfies the differential equation with the boundary 
condition: 



dp 



P(r) 



/(«o+ei+e 2 )-/(tio+ei) 


ei- 


/(«o+ei) 


-/(«o) 


e2-(ei+e 2 )Bi(p) 




/(«o+e2)-/(«o) 


ei — 


/(«0+eiH 


-ea)— /(«o+e2) 


e2 



= 1. 



t— >+oo 



This implies the following result. 



Theorem 1.5 TTie Cauchy problem for the equation (1.1), (1.5) in the sense of Definition L3 has a discontin- 
uous solution of the form 

2 

u*(x,t) =u a + ^2^e k H(- x + (j) k (t)) , 

k=l 



(1.31) 
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whe 



61 (t) = 0io (0) H t 

ei 

/(u + ei + e 2 ) - /(uo + e 2 ) /(u + e 2 ) - /(uo) 



»(*) 



ei e 2 

02o(O) H 1 

e 2 

f(uo + ei + e 2 ) - /(«o + e 2 ) /(u + e 2 ) - /(no) 



C2 



ei + e 2 



-(t-t*)H(t-t*), 



(1.32) 



ei 



e 2 



ei 



ei + e 2 



-(t-t*)H(t-t*) 



are phases. 

This solution of the initial value problem ( |1 . 1| ) , ( |1.5| ) is the unique generalized solution of the Cauchy probl 
in the sense of the integral identity (1.3). 



em 



This Theorem is proved in the same way as Corollary IA from Theorem LI and Theorem L3 
8. Description of the shock wave merging process. 

Thus, for t G (0, t*) two shock waves whose fronts move from the initial positions x\ = 0i(O), x% — 2 (O), 
where x\ < x\ along the straight lines 



with the velocities 



x\ = (j>w{t) 

X2 = 020 (*) 

d<fiw(t) 
dt 

dt 



x 1 /(ttQ+ei+e 2 )-/(« +e 2 ) , 
1 ' ei ' 

„0 I f (ua +e 2 )-f (up) . 
X 2 T e, I; 



/(Mo + ei+e 2 )-/(Mo + e 2 ) 
ei ' 

/(Mo+e 2 )-/(Mo) 
e 2 ' 



where 



d<j>io(t) ^ d<t>2o(t) 



dt ^ dt 

At the instant t — t* after the interaction shock waves merge constituting one new shock wave which for 
t G (t* , 00) propagates from the point (x* , t*) along the straight line 

x = x* + f{uo + ei+e2) - f{uo \ t-n 



ei + e 2 



at the velocity 



where 



dt 



■(*) 



f(u + e 1 + e 2 ) - f(u ) 



ei + e 2 



eie 2 



e 2 /(uo+ei+e 2 )-(ei+e 2 )/(uo+e 2 )+ei/(«o) 



/(«o+ei+e 2 )-/(uo+e 2 ) 



/(ti +e 2 )-/(M ) 



e 2 /(uo+ei+e 2 ) — (ei+e 2 )/(uo+e 2 )+ei/(tio) 



Remark 1.3 The same result for the Hopf equation can be obtained using the solution of the Burgers equation 



L B [u] = U t + (u 2 ) a 



eu. 



0, which is a regularization of the Hopf equation. It turns out that a simple 



solution of the Burgers equation can be found whose limit, as e — > +0, describes the process of merging two 
shock waves (see || 4.7.]). 

2 Proof of Lemma 1.1. 



First, consider the case when f(u) is an analytical function. We find the weak asymptotics of the product 

(2.33) 



j(x,a,s) =uj%i(— )uj% 2 ( X a ), 
e e 



which, according to Section 0, approximates the product of the distributions H m (—x)H n (—x + a), where 
m,n — 1,2,.... 
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It follows from the definition of the primitive of a distribution |^0|, Ch.I, §2.2.] that: 

J(a,e) = (g(x,a,E),<p(x)\ = -(g x (x, a, e), ^^(x)^, = J x (a, e) + J 2 (o, e), 
where (p(~^(x) = J"* 93(f) eif is the primitive of the function </?(x), 9?(x) 6 I?', J_ 00 99(f) (if = and 

Jfc(o,e) = -(^gk w (x,a,e),ip^(x)^, k = 1,2, 

/ \ 771 — 1/ ^ \ 77/ ^ "I" ^ \ 1 / ^* \ 

31a, (x, a, e) = -mw 01 ( )w 02 ( ), 

£ £ £ £ 

/ \ 777/ X . 1 X £1 . 

g 2 x{x,a,e) = -nu Q1 { )-u 2 { ). 

£ £ £ 

Making the change of variables x — eij we have 

Ji(a,£)=m J w™~ 1 {-t 1 )uj% 2 {-t 1 + -^{-^-^{e^dr,. 
Since the integrand decreases rapidly as — > 00, we have 

J 1 (a,e)=^ ( - 1) (0)m / ^(-ryK^-ry + ^^(-r?) dr? + 0( £ ), £ ^ +0. 
In an analogous way, changing variables x = a + er\ we derive the asymptotics 
J 2 (a,e) = ^3 2 (x,a,£),<^(x)^> 

= ^{a)n / ^(-^--HV^-^H^^ + O^), £^+0. 



Summing asymptotics (2.34) and ( p.35| ) we find the asymptotics 
J(a,e) = ^g(x,a,e),ip(x)^ 

= B™> n (^)<p(-V(0) + B^ n {-- s ) v { - 1 \a) + 0(e), e -> +0, 

where 

BT' n ( P ) = mfuj^i-rj^-rj + p^i-ridv, 
B 2 n ' n (-p) = nfu%(-T)-p)u% 2 l (-T))u2(-r))dr}. 

Here, as in the proof of Lemma 1.1, asymptotics ( 2.36| ) can be extended to the whole space D. 
Considering that 

^- x )(0)= f <p{§dti = {H{ r x),<p{x)), 



formula ( 2.36 ) can be rewritten in the weak sense as asymptotics ( 2.33 ): 



g(x, a, e) = B™< n (-)H(-x) + B™' n {--)H{-x + a) + O v ,(e), e - +0. 

£ £ 



(2.34) 



(2.35) 



(2.36) 
(2.37) 



(2.38) 



Using ( 2.33 ), ( B. 38| ), ( 1.20 ) and the binomial formula we find the weak asymptotics of the expression 



— X — X ~\~ CL "1 ™ 

d{x, t)«oi( ) + e 2 (x, t)u 02 { ) = e?(x, t)H(-x) + e"(x, t)H{-x + a) 

£ £ J 

+ B t (x, t, -)H(-x) + B 2 (x, t,--)H(-x + a) + Ot>>(s), £ -> +0. 

£ £ 



(2.39) 



Here, as it follows from ( [2.37 ), 

71-1 

B x {x,t,p) =Y,C*e k 1 (x,t)e%- k (x,t)B k 1 < n - k (p) 



k=l 
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fe=i 



n-l 



n— 1 



e\(x, ( — 77) drj. 



where are binomial coefficients. 

The last expression can be rewritten in the form 

Bi(x,t,p) = 

\ei(x,t)(jj i(-rf) + e 2 (x,t)uj 2(-ri + p)) - ( e t (x, t)u) i (rf) 
Analogously, we obtain 
B 2 (x,t,-p) = 

(ex(x,t)u 01 (-r) - p) + e 2 {x,t)uio2(-v)) - (^(x, t)uj 02 (-r])j 

Expanding the function f(u) into the Taylor series and using asymptotics ( 2.39| ) and the formulae for 
Bi(x, t, p), B 2 (x, t, —p), we fin d the weak asymptotics as e — ► +0 of f(uo + e\H\ {—x, e) + e2H2{—x + a, e))) , 
that is, obtain expansion ( |l.l7j ) and form ulae for the interaction switc h fun ctions ( 1.1 8[ ) . 

To prove the first formula from ( 1.1 9| ) we convert Bi(x, t, p) from ( |l,18|) to the form 



n-l 



e 2 (x,t)oj2(-T]) dr/- 



B 1 (x,t,p) = / f'(u + eiw i(-7?) + e 2 LU 2(~V + p)) eiWi(-rj) + e 2 w 2 (-?7 + p) 



dr) 

f(u + eiCJ i(-J?) + e 2 cj 2(-?7 + pj) - f'(u + e 2 cj 2(-?7 + p)) G2^2{-r\ + p) dr] 
- I f'(uo + eiU)oi{-rj))eiOJi(-r])dr} - / f'(u + e 2 uJ 02 (-V + p))e 2 w 2 (-77 + p) drj. 



Note that after the substitution — r\ + p — » —77' the second integral converses to the function B 2 (x,i, — p) 
from ( |1.18| ). Considering that lim r) ^ +oc (77) = 1, lim^-oo ujok(v) — we find that the first integral is equal 
to f(v,Q + ei + e 2 ) — /(iio) j while the third and the forth integrals are equal, respectively, to /(ito) — f{u$ + ei) 
and f(u ) - f(u Q + e 2 ). 

The first formula from (|l,19| ) is thus proved. 



Since lim p ^ +oc luq2{— r\ + p) — 1, relation ( 1.18 ) implies: 



Bx(x,t,oo)= lim B 1 (x,t,p)= f'(u + eiw Q1 {-rf) + e 2 ) - f'(u + e x u i(-r))) e\U\ (-77) drj. 

Integrating the last expression and taking into account that lim J( ^ +00 LOoi(rj) = 1, lim^-oo LOoiiv) — 0, we 
obtain 

B\ (x, t, 00)= lim Bi(x,t, p) = f(u + ei + e 2 ) - f(u + ei) - f(u + e 2 ) + f(uo), 

p— >+oo 



which is the second relation from (1.19). The other relations ( |1.19] ) are proved in an analogous way. 

The proof of Lemma is completed by using the Theorem on approximating a smooth function on a finite 
interval by analytic ones. 



3 Proof of Theorem 1.1. 



Substituting u'*(x, t) into the left hand side of equation (IT) and using ( 1.2C ) from Corollary we have 
Kt + [f( u *e)h = u ot + [f(uo)]x + e t H(~x + (j>) + cf> t e8(-x + (f) 
+ \f(u (x, t) + e(x, t)) - f(u (x, t))] H(-x + cP{t)) 
- f(u (x,t) + e(x,t)) - f(u Q (x,t))]5{-x + <f>{t)) + O v/ {s), e +0. 
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Thus, 



L[u*(x,t)] = L[u (x,t)] + ^L[u (x,t) + e(x,t)} - L[u (x, t)]\H (-x + <j)(t)) 
+ [<t>te{x, t) - [f(u (x, t) + e(x, t)) - f(u (x, *))] }<5(-x + <f>(t)) + O v >(s 



-0. 



(3.40) 



Setting the right hand side of ( |3.40| ) equal to zero, using Lemma [h^ (to separate the singularities) and 
Definitions 1.2, 1.3 of generalized solution, we derive Theorem 1.1 whose analog for strictly hyperbolic systems 
(when f(u) is a polynomial) was obtained in M, 8.1.]. 



4 Proof of Theorem 1.2. 



1. In the framework of our approach, to substitute the singular ansatz (1.15) into the Hopf equation means to 
substitute the smooth ansatz ( 1.16| ) u*(x, t, e) into the Hopf equation, and then to find a weak asymptotics, up 
to Oi)'(e), of the left h and side of the equation Lh [u*(x, t, e)] = 0. 
According to ( 1.21 ) , 



[H(-x,e)Y =H(-x)+Chy(e) 



-0. 



(4.41) 



To substitute ansatz ( 1.16 ) into the equation we use asymptotics ( 1.22 ) from Corollary L2 for the pair 
product H(—x, e)H(—x + a, e) a s e — > +0. 

In view of the aforesaid and ( 1.22 ) one should involve the dependence on e into shock wave phases, that is, 
seek them in the form <fi k = 4>k(t, s), k = 1,2. 

We cons ider the case when the shock waves amplitudes Ck{x,t) and the background uo(x,t) are constant. 
Using (pj~4i"l) - (|L23|) we find 



[u*(x, t, e)] 2 =ul+ (2u oei +e\ + 2e 1 e 2 B 1 (^^-)^H( -x + 4> x (t, e)) 



+ (2u e 2 + el + 2e ie2 B 2 ( - - x + 2 (i,e)) + Ov(e), e -» +0, 



where ip(t, e) 



^2 (i, — 0i(i,£) is the distance between shock wave fronts (discontinuities). 
Substituting the singular ansatz u*(x, t) and the asymptotics \u*(x, t)] 2 into the Hopf equation (1.2) we have 



L H [u*(x,t,e)] = 



fc=i 



e k —cj) k {t,E) - 2u e k -e%- 2 ei e 2 B k ({-l) k 1 p) 



S(-x + &k(t,e)) 



+ Ou'(e), e->+0, (4.42) 

where p = ^^' e ' is the independent variable of the interaction switch function and the estimate Ot» (e) in this 
representation is uniform with respect to the distance tp(t,e). 

From ( 4.42 ), using Lemma L2 (to separate the singularities), we derive the necessary and sufficient conditions 
for the relation Lh[u*(x, t, e)] = 0-d>(s) to be valid: 



dt 



<f> k (t,e) = 2u + e k + 2e 3 - k B k ((-l) k 1 p) 



1,2. 



(4.43) 



Before the interac tion of s hock waves, when i>(t,s) = (f> 2 (t,s) — </>i(t, e) > ce 1 "", where c > 0, < a < 1, 
according to Rem ark 1.1 and ( 1.24 ) we have, at least up to O(e), B\(^\ = 1 and B 2 { — = 0. Therefore, in 
this case system (4.42) becomes the system of equations which, according to Corollary 1.3 from Theorem [O], 
describes the dynamics of two non-interacting shock waves (taking into account that both the amplitudes e k 
and the background uq are constant): 



d(/)wt(t) 

dt, 
d(f> 2 ot{t) 

dt 



2u + ei + 2e 2 , 
2u + e 2 , 



(4.44) 



where 4>io(t), 4> 2 o(t) denote the phases of non-interacting shock waves. 
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Let ipo(t) — 4>2o(t) — 4>io(t) be the distance between the fronts of non- interacting shock waves. Define the 
interaction time t = t* as the solution of the equation ipo(t*) = 0. 

Equations ( [Oil ) describe the propagation of two shock waves before the interaction time, for t G [0, t*) 
and represent the Hugoniot conditions (1.27) along the discontinuity lines. Indeed, in our case for the lagging 
shock wave we have [u] = e%, [f(u)] = 2u§e\ + e\ + 2eie2, and for the advanced one we have [it] = e%, 
[f{u)\ = 2u Q e 2 + e%. 

Thus, before the interaction, i.e. for t G [0, t*), the shock waves move with velocities ( [1.44 ) along the 
straight lines 

si=<M*) = </>io(0) + (2w + ei + 2e 2 )t, 
X2= 020 (t) = ^20 (0) + (2u + e 2 )t, 

which intersect at the point with the following coordinates: 



(4.45) 



* = 02O(O)-01O(O) 

ei+e 2 ' 

* = 0io(O) + (2« o + ei + 2e 2 )t* 



where 0fco(O) are the initial positions of the discontinuities. 

2. To describe the interaction dynamics we shall seek the phases of shock waves as functions of the fast 
variable r = '^ } t ' G M and the slow one t > 0: 



4>k{t,e) d = f (p k (r,t) = fe(t)+^oWfe(T), 



(4.46) 



where the functions (f>ko(t) defined by equations (4.45), for t G [0, t*), are extended by these equations (4.45) 
for all t € [t* , +oo). If t > then t < t*, i.e. the interaction has not occurred yet; if r < then t > t* , i.e. the 
interaction has occurred. 

Since we consider piecewise constant solutions, we can take the velocities of all the shock waves formed 
before and after the interaction to be constant, and seek the perturbations for the phases 4>ko{t) as functions 
0fci (t) only depending on r. 

We set the following boundary values: 



4>ki (t) 
#fei(r) 



T— >— OO 



0. 



(4.47) 



That is, the derivatives of the phases with respect to the fast variable r tend to zero as |t| — > oo, while the 
phases themselves tend to zero as r — * oo. 

Finding the limit values of the perturbations as r — > — oo 



5fcl. 



we find the phase limit values 



4>k{r,t) 



= 4> k ,-(t) = (j>ko{t) +4>a{t)<t>kt,- 



and thus define " the result" of the interaction of shock waves for t > t* 
Let ipi(r) = 021 (t) — 0n ( T ), then the full phase difference is 

^(t,e)=^(t)(l + Vi(T)), 



the independent variable of the interaction switch function ( 1.25 ) has the form: 

p(t) = ^M=t(1 + V>i(t)). 



(4.48) 



The phase derivatives with respect to time are represented by the following relations: 
d a a r\ def d 2 ( A d Mt) d^ (t) d 



(4.49) 
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Taking into account the boundary conditions (4.47), we obtain the limit values of the phases and their 
derivatives with respect to time asr-> — oo: 



4>k,-(t) = <pko(t) + lpo(t)(j)kl-, 



(4.50) 



Substituting (4.46), (4.44) into ( 4.43| ), we derive, for all t > and r 6 (— oo, oo), the following system of 
equations with the boundary conditions ( |4.47| ): 



d<p w (t) dip (t) d 



dt dt dr 

d<fao(t) di>o(t) d 



dt 



dt dr 



T<f>n(r) 
T(j>2i (r) 



2u + e 1 + 2e 2 B 1 (p), 
2u + e 2 + 2e 1 B 2 {-p), 



(4.51) 



where <pko{t) are the functions defined, for all t >,0 by equations (4.45). 

Subtracting one of the equations (4.51) from the other we find the differential equation satisfied by the 
function p(r): 

Pr = F{p), 



T — ^ + 00 



= 1. 



(4.52) 



where the boundary condition ( 4. 52 ) follows from ( 4.47 ) . 

Taking into account the first relation (1.24) and the fact that, according to (4.44), V'ot(i) = — ^ei + e 2 ^j, we 
find the right hand side of the differential equation 



HP) 



e 2 (l-2Si(p)) - ei (l-2B 2 (-p))] = 2B x {p)-l. 



(4.53) 



According to Corollary |l.2| , the function -Bi(p) is smooth and its values belong to the interval (0, 1). 
Therefore, independent of the choice of the approximating functions Wfc(£), k = 1, 2, the equation F(p) = has 
at least one root po which solves the equation B\(p) = 1/2. 

The differential equation (4.52), (4.53) is autonomous. For this equation the following statement holds. 



Proposition 4.1 For the autonomous equation 



with the smooth right hand side F(p) to have a solution such that 



Pill 

T 

P(r) 



T—^-\-OC 



1, 

Po, 



(4.54) 



(4.55) 



where po is a constant, it is necessary and sufficient that the following conditions hold: 



F(p) 



= 1, 

F(p ) = 0, 
F(p) > 0forp>p , 



(4.56) 



where po is the maximal root of the equation F{p) = 0. 

In addition, if po is an ordinary (nonmultiple) root of the equation F(p) = then p—po = 0(t~ n ), t — > — oo 
for all N = 1,2,.... 



Proof. Suppose that the limit relations ( 4.55|) hold, and f ,f + 1 € (po, +oo). Then, integrating ( 4.54 ), we 



have 



p(f+l)-p(f) 



F{p{r + f )) dr. 
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Since according to the second relation (4.53), the left hand side of this equality tends to zero as f 
limit of its right hand side is also equal to zero: 



the 



lim 



F{p(r + f )) dr 



f F(p )dr = F(p ) = 0. 
Jq 



That is, po is a root of the equation F(pp ) = 0. 

It follows from the first relation ( 4.53| ) that F(p) > for p > po, and consequently, po is the maximal root. 
Since for n, r S (pa, +oo) 

p(t) _ pin) l 



F(p(r'))dr', 



(4.57) 



passing in (4.57) to the limit as r - 
rule, we prove tha t lim ^ +oo F(p) 



-co, taking into account the first relation (4.55) and using the L'Hospital 



Conversely, if ( 4. 56] ) holds then p = po is a solution of equation ( 4.54 ). Integrating ( 4.54 ), we have 

rP dp' 



F(p'Y 



Since the function F(p) is smooth, the function l/F(p) has a non-integrable singularity at the point p = p , 
and the integral jFpp\ diverges as p — > p . 

That is, as p approaches po, we have r — > — oo and therefore, p — po is an asymptote for all integral curves 
within the region po < p < +oo. 

The first relation ( 4.55 ) follows from ( 4.57 ). If po is a simple root of the equation F(p) — then F(p) — 
(P - Po)g(p)), where g(p) > for p > p . 

The second statement is proved by passing from the differential equation ( 4.54 ) to the differential inequality 

< (p — po)M, where M = max p > Po g(p). Integrating this inequality, we find that p — po < Ke T . Thus, 
p — * po more rapidly than any power of | -r | 1 as r— > — oo. 

The proposition is proved. 

— oo, we have p = r(l + ?/>i(t)) 



By Proposition 4.1 



as r 



po, here Bi(p) -> Bi(p 



1/2 and 



B 2 (-p) - B 2 (-p ) = 1/2. 

Thus, passing to the limit in (4.51), as r — + — oo, and taking into account that, according to (4.47), the phase 
derivatives with respect to the fast variable r tend to zero, we derive the following limit system of equations 
which describes the shock wave evolution after the interaction, i.e. for t > t*: 



d(j>io(t) dip (t) 

dt dt ' 

d(/)2o{t) , dipo{t) 



dt 



dt 



ii.- = 2u + e 1 + e 2 , 
2i.- = 2u + e 2 + ei. 



(4.58) 



Remind that here the functions 4>ko(t) ar e determined by equations ( 4.45| ) for alH > 0. 
It is clear from (4.58) that phase limit values of both shock waves coincide, that is, 



,_(*) = &,_(*) W<Mt). 



(4.59) 



Thus, it follows from (|4.58|) , ( |4.59|) that after the interaction, i.e. for t > t*, the discontinuities merge 
together constituting a new shock wave whose dynamics is determined by the equation 



d4>-(t) 
df~ 



= 2wn 



ei + e 2 



(4.60) 



with the initial value 0-(t*) = <f>ko{t*)- 

Equation (4.60) represents the Hugoniot condition (1.27) along the discontinuity line after the interaction, 
since [u] = Uo — (u + ex + e 2 ) = -(ej + e 2 ) and [f(u)j = (u ) 2 - (u + e x + e 2 ) 2 . 

3. Now we describe the dynamics of the shock wave merging process. To this end, substituting the phases 
<t>ko(t) from (4.44) into (4.51), and using the relations B\(p) + -B 2 (— p) = 1 from (1.24) and ipot(t) — — ^ei + e 2 ^, 
we derive the equations determining the perturbations for the phases </>fei(r): 

, fc _i 2e 3 _ fc 



d_ 



T^ii(r) 



(~1) A 



ei + e 2 



(l-Bi(p)), fc = l,2. 
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Integrating these equations, we have the following expression for the phase perturbations 



fcito = (-l) 



k-l 



2e 3 _ fe 



l-B 1 (p(T , )))dT', k = l,2, 



(ei + e 2 )r J 

where p — p(r) solves the differential equation with the boundary condition ( 4. 52] ), (4.53): 

dp 

P(r) 



(4.61) 



dr 



= 2fli(p)-l, 



T — >+00 



Verify that the functions <^fei(r) found in (4.61) satisfy the required properties. 

We have chosen the approximating functions uik(z) so that they either have compact supports or decrease 
sufficiently rapidly as \z\ — > oo, for example, |wfc(z)| < Cfc(l + |z|) -iV , N — 1,2, . . ., fc = 1,2. Therefore, we 



have (a>i *W2)(?j) < K(l + \z\)~ 
p^+oo, Bi(p) = <D(p~ N ) as p 



as 2 



JV = 1,2, 



and B x {p) = (a>! * wa)^) rfr? = 1 + 0(p- JV ) 



as 



If t — > oo then p — > oo and, according to (4.61) and Remark 1.1, we have 



4>ki{r) 
dr 



0(0, 
0(t-1). 



Note that p = po is an ordinary (nonmultiple) root of the right -hand side of the differential equation ( 4.52 ), 



(4.53). Indeed, F(p ) = 2B 1 {p ) - 1 = and accordingly to (1.23), 



F '(Po) = 2 (^i * w 2 )(p ) 



Wfc(0^3-fc0o + d C > o 



since > 0, fc = 1, 2. 



Thus, if r — > — oo then p — ► po and according to Proposition 4.1, p — po = 0(\t\ ) for all 2V = 1,2, 



Therefore, using Remark 0, we have Bx(p) = 1/2 + 0(|r| 



-oo. According to (4.61), we obtain 



&i(r) = (-l) fe - 1 ^ + 0(r- 1 ), 
V-i(r) = -l + O^- 1 ). 

Calculate the limit of the expression t HmSi1 as T _> — oo. To this end, rewrite this expression in the form 



r #«(r) =( _ 1)fe -x_f3-^_ (l _ jBi(p)) 



rfr ' ei + e 2 

Applying the L'Hospital rule to the limit in square brackets we find 

lim T(f) klT (T) = 0. 



r Q {i-B l { P (T')))dT> 

1 r(l-Bi(p)) 



As one can see, the functions (r) obtained in (4.61) have all the properties required. 
Thus, we have proved the Theorem 1.2. 



5 Proof of Theorem 1.4. 



1. We consider the case when uo(x,t) an d efc (x, t) in ansatz ( 1.15| ) are constants. In this case, therefore the 
interaction switch fun ctions from Lemma IT Bk(x,t,p) = Bk(p). As in Section |], we shall seek shock wave 
phases in ansatz ( |l.l5| ) as functions of e, i.e. 4>k = s)> k = 1,2. Denote the distance between shock wave 
fronts by (t/j(t,e) = (f>z{t,e) — <fii(t,e) 



Substituting the singular ansatz (1.15) and asymptotics (1.17), given by Lemma 
have 

L[u*(x,t,e)]=u* et + [/«)], 



1.1 



into equation (1.1) we 
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fe=i 



efe- 



dt 



/(«o + e fc ) - /(m )] - S fc ((-l) fc - V) + e)) 



+ Od/(e), e->+0. (5.62) 

where p = ^(*' e ) and the estimate Ou'(e) in this representation is uniform with respect to the distance ip(t, e). 

From ( p% using Lemma [T7^ , we derive the necessary and sufficient conditions for the relation L[u*(x, <)] = 
Ot>> (e) to be valid: 

/(tt + ei) - /(u ) 



ei—(l)i(t,e) = 
at 
a 

e 2 — 2 (i,£) = 



f{u + e 2 ) - f(u ) 



Si(p), 
B 2 (-p). 



(5.63) 



When 



V>(*,e) = 2 (t,e) - <M*> £ ) > ce 



where c > 0, < a < 1, that is, before the interaction of shock waves, according to Remark 1.1 and (1.19), we 
have, up to 0(e): 



*(£) 



f(u + ei + e 2 ) - /(ito + ei) - f(u Q + e 2 ) + /(m ), 



B 2 (--) = 0. 



rem 



It follows from the above that system (5.63) turns into the system of equations which, according to Thco- 

iteracting shock waves: 

= f(uo + ei + e 2 ) - /(ito + e 2 ), 



1.1, describes the dynamics of two noninteracting shock waves: 

io(t) 



ei 



e 2 - 



dt 
d (f>2o(t) 

dt 



f(u + e 2 ) - /(ito). 



(5.64) 



Equations ( 5.64 ) obtained above represent the Hugoniot conditions ( 1.27 ) along the lines: [u] = e\, [/(it)] = 
f(v,Q + e\ + e 2 ) — f(uo + e 2 ) for the lagging shock wave, and [u] = e 2 , [/(it)] = /(ito + e 2 ) — /(ito) for the 
advanced one. 

Denote by ipo(t) = 4>2o(t) — 4>w(t) the distance between the shock waves before the instant of the interaction 
t = t* which we define as a solution of the equation ipo(t*) = 0. 



Thus, before the instant of the interaction, as t e (0, t*), the shock waves move with velocities (5.64) along 
the straight lines 



XI = 010 (*) = X\ + nuo+ei+e 2 J- n u +e 2 ) ^ 

which intersect at a point with the following coordinates: 



(5.65) 



f 



eie 2 



e2/(«o+ei+e 2 )-(ei+e 2 )/(tio+e2)+ei/(«o) ' 



f(u +e 1 +e 2 )-f(u +e2) 



e 2 x 2 - 



f(u +e 2 )-f(u ) 



e2/(«o+ei+e2)-(ei+e 2 )/(Mo+e2)+ei/(Mo) 



(5.66) 



where (f>ko(0) are the initial positions of the discontinuities. 

2. As in Section |^, define shock wave phases by formula ( 4.46] ) where the functions 4>ko(t) determined 
by equations (5.65), as t S (0, i*), are expanded by the same equations ( 5.65| ) for t G (t* , +oo), as was done 
previously, and the phase derivatives with respect to time are given by formulae ( |4.4£ ) . The full phase difference 
is 

1 p(t,s)=Mt){l + Mr)), 



where ipi( T ) — < / ) 2i(' r ) — <Pii( T )- Here the indepen dent variable of the interaction switch function is set in ( 4.48 ). 

For <f>ki( T ) we se t the boundary conditions ( 4.47 ). The limit values of phases and their derivatives with 
respect to time, as r — * — oo, were determined in (4.5C). 
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Substituting ( 4.46| ), (4.49) into ( |5.63| ), we obtain, for all t > and r e (—00, 00), the system of equations 
with boundary conditions (4.47): 



ei 



dt dt dr 

(d<j) 20 {t) dipo(t) d 
e 2 ■ 



di 



dt dr 



t<Pu(t) 
T<p2i (r) 



= f(u + ei)-f(u ) + B 1 (p), 
) = f{u + e 2 )-f(u ) + B 2 (-p), 



(5.67) 



where 4>ko(t) are the extensions of the functions 4>ko(t ) determined by equations ( 5.64| ) for all t > t*, the 
extension being determined by the same equations (5.64). 

Subtracting one of the equations in system ( |5.67| ) from the other we obtain the following differential equation 
with the boundary condition for p(r): 

Pr = F(fi), 

P 



T — >' + 00 



= 1. 



where 



Hp) 



dt 



f(u + ea) - /(«o) + B 2 (-p) f{u + e x ) - f(u ) + B x {p) 



C-2 



ei 



(5.68) 
(5.69) 



here the bou ndary condition follows from ( [4.47| ). 
Formula ( 5.64 ) implies: 



tpot(t) = 



f(uo + e 2 ) - f(u ) f(u + ei + e 2 ) - f{u + e 2 ) 



('2 



ei 



(5.70) 



Taking into account the first correlation from (1.19), convert function (5.69) to the form 

Hp) 



f(u + ei+ e 2 ) - f(u + ei) 


ei - 


f(u + ex) - f(u ) 


e 2 - (ei + e 2 )Bi(p) 




f(u Q + e 2 ) - f(u ) 


ei - 


f(u + ei + e 2 ) - f{u + e 2 ) 


e 2 



(5.71) 



Using th e lim it values of the switch function _Bi(±oo) from Lemma |1.1| we find the limit values of the function 
F{p) from ( |o~7ll ): 



F{P) 
F{p) 

Rewrite the limit value F(p) 
F(P) 



= 1, 



/(ti +ei+e 2 )-/(tio+ei) 


ei- 


/(«o+ei)-/(M ) 




/(n +e 2 )-/(«o) 


ei- 


/( u +ei +e 2 )-f (u +e 2 ) 


e2 



(5.72) 



p — >— 00 



from (5.72) in the form 



f(u a + ei + e 2 )ei - f(u + ei)(ei + e 2 ) + f(u a )e 2 
f{u Q + ei + e 2 )e 2 - f(u + e 2 )(ei + e 2 ) + f{u )ei ' 



In view of our assumptions, the function f{u) is convex downwards, that is, f"(u) > on the range of 
the solution it, and the amplitudes are positive, e%, e 2 > 0. Under these assumptions, both the numerator 
and denominator of the last fraction are strictly positive, since for a function convex downwards the inequality 
(x 2 — x)f{x\) + (xi — x 2 )f(x) + {x — a^i)/(a; 2 ) > holds for x\ < x < x 2 . In our case, x\ — uq, x = uq + e\, 
x 2 =u a + ei+e 2 . 

Thus, as p — > ±00, the limit values of the right hand side F(p) of the differential equation (5.68) have 
opposite signs. 

According to (1.18), the derivative of the interaction switch function has the form 
d 



-j^Fi(p) = e i e 2 / f"(u + eiujQi(-r)) + e 2 u 02 (-r) + p))ui{-rj)w 2 {-r) + p) dr\. 



It is positive since f"{u) > 0, Wfc(?7) > and e±, e 2 > 0. It follows from here and from (5.71) that the derivative 



——F(p) = ( e i + e 2) 

dp f{u + ei + e 2 )e 2 - f(u + e 2 )(ei + e 2 ) + /(w )ei dp 



Bi{p) 



(5.73) 
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is positive, and therefore F(p) is an increasing function. 

Since the function F(p) is smooth, the equation F{p) = has a root pq. Then, according to Proposition 4.1, 
the solution p of equation ( 5.68) ) tends to po, as r — > — oo. 

Due to the smoothness of the function F(p) and the aforesaid, the equation F(p) = has at least one root 
Pq. Let pq be the maximal root of the equation F(p) = 0. Then, according to Proposition 4.1, as r — > — oo, the 
solution of equation Q5.68 ), p tends to pq. 

Pass to the limit in ( [5.67 ), as r — > — oo, taking into account that by (4.47), r0fci r (r) — > 0. From here we 
derive a system of equation describing the evolution of shock waves after the interaction, for t > t*\ 



d<j>ig(t) 

dt 
d(f>2o(t) 

dt 



dlppjt) _ f(u +e 1 )-f(u )+B 1 (po) 

dt 61 ' 

dj>o(t) _ f(u +e 2 )-f(u )+B 2 (-p a ) 

dt ~ 62 



As shown above, we have p = rfl + ^i(t)) — ► po as r — ► — oo. It follows that ^i(t) 
V-i(t) = -l + 0(r- 1 ). 

Since s/>o(t) = <p2o(t) — </>io(i) and f/'o(i)0(r _:L ) = 0(e), we have asr^ — oo: 

Mr) = 0ao(t)+V^)(*)0ai(r) 
= 02o(i) + M*) f^iiW - 1 + Oir- 1 )) = Mr) + 0(e). 



(5.74) 



— 1, that is 



That is, up to 0(e) we have 



def 



02,-(*) = &,-(*) = 0-(t). 



(5.75) 



It is clear from (5.74), (5.75) that after the interaction, i.e. for i > t*, discontinuities merge constituting a 
new shock wave whose dynamics is determined by the equation 



d$-(t) _ f(u +ei) - /(u ) + Bi(po) 



ei 



(5.76) 



Here the following relation holds 

f(u + ei) - /(u ) +Si(p ) /(u + e 2 ) - /(ito) 



S 2 



Po) 



ei 



f'2 



(5.77) 



From ( 5.77 ) using the first relation from ( 1.19 ), by Lemma 1.1 , we find the value 



Bi(po) = 



f(u + a + e 2 ) - f(u + ei) ei - /(u + ei) - f(u Q ) 



e2 



ei + e 2 



(5.78) 



Substituting value Q5.78 ) into (5.7£) we obtain the equation which describes the propagation of the shock 
wave resulting from merging the initial waves: 



• (*) _ /(«o + ei +e 2 ) -/(uo) 



ei + e 2 



(5.79) 



with the initial value <fr-(t*) = <f>ko(t*) — x * ■ 

Equation (5.79) is the Hugoniot condition ( |1 ,27| ) along the line of discontinuity after merging the shock 
waves, since [u] = u - (u Q + e\ + e 2 ) = — (ei + e 2 ) and [f(u)} = f(u + e\ + e 2 ) - f(uo)- 

3. Consider the dynamics of the proces s of m erging two shock waves into a new one. Substitute the p hase s 
of noninteracting shock waves 4>ko (t) from (HH) into (|j37]). Fro m here, using the first correlation from ( |l,19| ) 
for the interaction switch functions Bk(p) and expression ( 5.64 ) for ipot(t) we obtain equations to determine 
0*1 (r), k = 1,2: 



d 

d7 



T^ki(r) 



24 



/(m + ei + e 2 ) - /(u + ei) - /(u + e 2 ) + /(«o) - Bi(p(r)) 



C2 



/(«0 + e 2) - /(uo) 



/(u + e! + e 2 ) - f(u + e 2 ) 
Integrating equations (5.80) we derive the following expression for 4>ki(j): 
4>ki(r) 

" /(u + ei + e 2 ) - f(u + ei ) - f(u + e 2 ) + f(u ) - Biipir')) 



(5.80) 



ei 



("If" 1 



e3-fc 



/(u + ei + e 2 ) - f(u + e 2 ) 



e 2 



/(u + e 2 ) - /(u ) 



dr' 



(5.81) 



ei 



where p = p(r) solves the differential equation with boundary condition ( 5.68| ), fl5.69| ). 

As in the case of the Hopf equation, it can be verified that the functions 0fci(r) found in ( 5.81 ) satisfy the 
presupposed properties. 

Let r — * oo, then p — ► oo, p ~ r and, according to ( 5.81 ) and Remark 



0ki (r) 
^fci(r) 



0(0, 
0(0. 



According to (5.78), p — po is an ordinary (nonmultiple) root of the right-hand side of the differential 
equation ( |5.68 ), ( 5.71 ), i.e. F(p ) = 0. But in view of to ( 5.73 ), j^F(p) > 0. Consequently, p = p is an 
ordinary (nonmultiple) root of the right-hand side of the differential equation ( 5.68 ), ( 5.71 ). 

Thus, if r — * — oo, then p — ► po and B\(p) tends to the limit value Bi(po) given by formula (5.78). 

p — po = 0(\t\~ n ) for all N = 1,2,.... Thus, using Remark [El], we have 
— > — oo. Therefore, ( 5.80| ), ( ^.81 ) imply the same estimates as in the case of the 



According to Proposition 4.1 
B 1 (p)=B 1 (p ) + O(\r\- N ), ■ 
Hopf equation: 



fcl (r) = (-l) fe -^ + 0{t-% 
Mr) = -l + O^" 1 ), 
_d4>ki(r) 



dr 



0(t-1) 



Thus, the functions <pki( T ) constructed in (5.81) have all presupposed properties. 
Thus, we have proved the Theorem 1.4. 

The authors are greatly indebted to E. Yu. Panov for many fruitful discussions. 
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